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Forecasts of ARIMA processes are generally 
Equation form. This is the approach favoured by 
subsequent authors. The purpose of this note is 
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In their bcck ; Box and Jenkins (1970) strong -7 recommend tuat forecasts 
of ARIMA processes be once using the Difference Equation form because it is 
the simplest approach. All subsequent textbooks have endorsed tnis view to 
the extent that very few of them consider forecasting in the Integrated Bern 
in any detail, A notable exception to this is the recent book by Abraham an 
Ledolter (1983) which contains seme useful detailed discussion of one role c 
the Eventual Forecast Function (EFF) in generating forecasts. The purpose of 
this note is to indicate that the Integrated Form does have value and is v/or 
of consideration. 7/e begin with a brief review of the derivation of she 
necessary components for this form of the forecast. 

Suppose {X + } is an ARX1A process satisfying a!3)X. = 

u ~ t 

where B is the Backshift bee ra tor whicn is defined bv 3~X\, = 



= O i ~ 
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:ed in the c - ana 6 - or era tors. For lead time 



l cu.Xj.lT-k) = D , where X^i) is the forecast mace at time t 

K - 0 

X . and takes the value X . if i < 0 . The solution cf this :iffer c 
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m oe written as 

o 

Y t rp \ _ O , v p / m • 



where f, (T) are deterministic functions of T and may include p o Igno on 
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exponentials, sinusoids and products cf these. They nay also ce dunrpy 
variables generating a seasonal pattern. 

As a r^ ore sen tat 1 on cf the T— steo ahea i forecast vl ' a cove is vails 
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./here {e is the sequence of one-step .ahead forecast errors, i.e. 
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Box and Jenkins (1970) show that the current value of b^ = , . . .b ) 

may be obtained from To via a linear equation b = Lb + he . The 

L* -i- T> O -i- o 



matrix L effects the changes in the coefficients in revising the time origin 
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from { t — 1 ) to t , and can be obtained as L = F J( F,, +1 , where F, f is the (p*p) 
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matrix v/i th ( i , j ' oh element > The vector h = , v/h 
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and b, is the coefficient of 
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;(B) = £ ( B ) / -x '» B J , tr.e usual moving-average representation of the ARI.MA 



trocess. 

V/e ray note that the revision matrix L is block diagonal. Each block 



corresponds to a real (repealed) linear factor or conjugate pairs of complex 



factors in a(B) . The effect of this structure is that individual components 
can be ...cr.it ore d, revised and projected independently of the others. Thus, for 
example, any linear tr-euu ar. i seasonal factors can be obtained at each time t 



given only their values at time (t-1) and the latest forecast error e A . 

o 

Details of the derivation of the revision equations for the linear trend and 
seasonally differenced Ah I MA model forecasts are given in McKenzie (1934b). 

They ran be obtained .vithcut direct evaluation and inversion of the matrix F . 

If ve /is:, to genera t-e r.at merely individual components of the SFF but the 
forecast itself using * r.e Integrated for:... v/e must use (1) or (2) above. If 
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The derivation of this result is given in an Appendix at the r- •. • , ... 

There it is shown that d„, . = d . , (j = 1,2,...,\I-T ; T = 1,2,...,M). 

Thus, it is necessary to obtain only {a^ ^ : T = 1,2,...,'!}, and these are 

easily derived from the following .matrix -.juation. 
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The triangular nature of these cpn: out evidently takes the.'., particularly ea 
to solve. This is especially true when ,ve note that ’ ! = q - p is very rare 



large . 
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Bex and Jenkins give sene non -seasonal examples in tneir 
scnal one is giver, in ' k'henz. ie !l?:en! e repressing the .veil -known airline .to da 
forecasts in c:rr=nt level, yndi -;r.t and seasonal factor form. V.e give one 
other brief example here * j illustrate the evaluation of the d-coef f icients . 
Consider 

(1 - rial - BiXj. = (1 - 8&‘ T )a 4 . . 
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The Integrated Terra of the forecast is now 
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(ii) Forecast Equation 



x.(k) =o + kb; 
t o 1 



k > 1 



(5) 



There is a fixed computational investment in generating t = 



t v t 
'o'~I 



before any forecast can be obtained. Once this is achieved, ho.vever, 
forecasts for any lead times are easily derived using equation (5;. All tne 
storage requirements relate to the revision equations (A), and here involve 

\-l’ V h 2 ar ‘ d V 

For the Difference Equation form, the forecasts are generated directly: 



X, ( 1 ) = X. , 2) + *,• 
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Note that there are no revision aquatic 



net venerate forecasts. 



However, it is clear that tne forecasts of lead fines h-1,2 and 3 play a role 
in equations (o) and (7, 1 similar to that of iy in equations and (5!. 

I 

Another important point nere is that to obtain forecasts for any lead time 
using (6) and (7j we must first generate the forecasts for all shorter lead 
times . 

To illustrate the differences between the two approaches we consider some 
simple forecasting scenarios. Suppose wo wish to forecast for lead times 1: = 
1,2,...,T. The computational requirements are clearly comparable. Tne 
Integrates Form requires less storage to generate the forecasts, but the 
revision -equations necessitate a little more arithmetic tnan the use of (c) 
and (7). Eurtcse now tnat the leas times of interest are not conS'.o stive . 



For example, sunt 



s mg one mt- 



ts for lead times h- 1,2 , 3,0 



b has been obtained. However, the Difference Equation form requires the 
generation of ail 13 forecasts. Alternately, we may routinely forecast for 
lead times 1,2 and 2 and occasionally require forecasts for others, e.g. 6,12, 
13 etc. Again, it is easier to generate these via (*+),(?) than (6), (7). 

Another prediction cf common interest is the cumulative forecast, i.e. 



the forecast of Y_CO = 



1 X (k). Using the 



l X . It is given by Y O'! 

Difference Equation form we must generate X. (k) for k=l,2,...,T and then sun 
them. For the Integrated Form, we can use a single forecast:- 
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. r i e ty c: • ore cas tint croDler ?.s are 



thus rcre efficiently treated by this approach* To assess the cost cf such an 
investment , note that the tir.er.sion f the vector b^ is the order of the 
autcre tressive- :if Terence operator in the AT I '.'A process, e . g . (p+d) in an 
ARI.V.A !p , d , : ; process. Tnus, revision ;f is rcugnly comparable to the 
o rouu 1 1 i tn of forecasts for nae ’ i^st o^t ^ l^^d tones i i oi r ' 1 ’ the p f *' p x*e n c e 
Ecuati :n fern. Tlearly, if no lead tines beyond (p+d) are to be predicted, 
the liffer ence creation for:?, is ncr a eff ioie.nt in general. However, when, lead 
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r^aca j^C'^rus , as v/e hate 



seen, n 



On toe other hand, there is a most Important advantage en 0 yj^ l.// the 
Integrated form. It lies in the area cf interpretation. This is an aspect s:' 
forecasting which should not be underestimated. V/e may recall Item's 
comments (1974) that managers who require forecasts may be prepared to accept 
trends and seasonal effects because these correspond to familiar ideas. This 
view is merely one of many such cr is de cceur still heard by forecasters. The 
true value of a forecast is not invested solely in its accuracy but a is o in 
its credibility. The former can be assessed only after the event forecast, 
whereas the latter will determine whether the forecast is used at all. TL*= 
modelling of tine -series as AKIM A processes is now a common practice, thanks 
to the availability of a variety :f powerful computer packages. fevers: el - 
a major hurdle for most osers is still tne practical interpretation cf the 
f orecasts . 

The reason that the Integrated form enjoys an advantage here is that in 
generates the forecast via the Eventual Forecast Function IFF . This is a 
linear combination cf deterministic functions of the lead time T. For fi;n : 

T, the functions are sum exactly, but their coefficients in the li:oar 
combination adapt v/it.n each new observation. Thus , in the AEI'iA' i ,2 , 1 ) 
example, the EFF is given by equation !;). The deterministic functions are 
f 1 (T)=l and T)=T. Their coef f icients , the components of b^, are revised 
with each new observation via equations fi.'. The new observations allow the 
forecast to adapt, while preserving the basic structure of the EFF. 

This structure is most important because the individual deterministic 
f -.unctions union appear in tne EFF car. often be readily interpreted to toms 
of familiar concerts such is trams, srrwta an i seasonality. m is 
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jus. suer: a linear combination of deterministic functions. Surely , the 
success and '.vide acceptability of these systems is at least partly due to the 
fact that their predictions are cast in a familiar form. In addition, more 
recently popularised systems such as the Bayesian models of Harrison and 
Stephens (19 7 6), cr the Kalman Filtering approach discussed by Harvey (1983), 
are also constructed around these same familiar com arenas of trend and 



seasonality, etc. Cn the other 
modelling and forecasting ARIMA 



hand, none of the software available fc 



:rc cesses can generate forecasts in this 



readily understood form. Nevertheless, such interpretations are available via 
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revision equations yield current values for the non-seasonal corn .rent, tq 
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(b ,b. ), and the seasonal component 3 X = (3, ,S~ , . . . ,3, ^ ) 
o ' 1 * — t x 2 ' 9 12 

have useful and readily understood interpretations* In the non -seasonal 



component, b is the current level of the 
7 o 



>eess, and b, the current gradient 



of the linear trend, i.e. the predicted rate of change in 

the seasonal component, is the additive seasonal factor for rrcr.th (t+I:). 
predicts the amount by which the data will deviate from the process level h 
months from new 'time t). 
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quote the non-seascr.al component without reference to the seasonal one but,, 
from equations (S) and (9), we can calculate them separately. Thus, for 
example, we could forecast using the Difference Equation form and monitor the 
level and trend directly via (3). 

Exhibit 2 about here 

The 'values of the seasonal component are revised with every observation 

same way as the ncn-seasonal one. Since it involves 
, in fact, very stable, they are not reproduced here. 
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and can be plotted in the 
twelve maths and they are 



It is just as revealing, in thi. 



- b‘. 

c 



is dismlayed in Exhibit 



n rr >y 






• a ror tr S 



manUSI 



always high and approximately the same. The ether months are equally 

‘xc-mtion of lecember, which appears to ce on a 



U -a -a, * o tr 



to not • 



-r 1 " - 



and fioo 



1 070 > 












) those authors 



, 1 0 y 



of this forecast 



j ^ r 4 o ■> 



sent i : , 



it it nas sirse ;ter 



n ventun x “ c reca s t n 



ctMMM 

garer is to raake a case for the usefulness of the 
Mm forecast. A brief rev lev/ of its derivation is 

tart of it is given in 

iv ^ h ^ * q t ’ r * *r> .v 1 i + p r»o - 1 r»p> 

rom its direct use of ' 






when a subset of these data ( ur 



•' m o o •*- 






nave erroo.asxs 



iced two coints. 






First, it need not be, as is often supposed, computations lip inferior 
the Difference Equation approach. Secondly, it offers, via the IF l , 
considerable advantages for interpretation of the forecast. In both cases, 
argument is illustrated with an example. 
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Derivation of the d-coefficlents 

First, note that substituting equation (2) into the usual forecast 



revision identity 

- ^.- + l e t 



X + (T) = X + _, ( T+l ) + i 



yields d„ = d„ 



'h (-* 

> 



i _1 



it is necessary to obtain only (cL 



and using the revision equation 

( i =1 9 m = ] O \i\ 

'J x y - 1 ; *> ; • • • ' ' 

: T = 1,2,... ,1!} and this .nay be 



Applying the usual conditional expectation arguments to X. to derive T-step 

t 



i v- oq a 



ists : or : = 



■ y • • • ; i > 



yields 



a- .o-r-.-n 



f ^ ~ ^ o vr 



ecus 



a, : . ; . 



- r; o ' O * 



-nans sice ot 



nay ce written 



x < 



-a ) i 



v/ne re 



t urther 



eq union yie.oc 



11 ) 






.ve i m 



, r r + 1 



■ v ^ ; .r 



< i < 



1 < " < .7 



:n 



- . i e ; id . 



r: uati :n 



v^> 



REFERENCES 

Abraham, E. and Ledolter, J. Statistical Methods for Forecasting, hu.v York 
John Wiley, 1983. 



Box, G. E. P. and Jenkins, G. M . Tine series analysis: forecasting and 

control , San Francisco, Holden-Day, 197C. 

Godolphin, E. J. "A direct basic form for predictors of autcre 
integrated moving average processes", Bio.ne trike , 62, 2, 1 

Harrison, P. J. and Stephens, C. F. "Bayesian Forecasting", J. R. Statist. 
Soc. B , 33, 3, 205-247, 1976. 

Harvey, A. C. "A unified view cf statistical forecasting procedures", 

J. of Forecasting , to appear 193k. 

Hillner, S.C. and Tiac, G.C. "Likelihood Function cf Stationary 
Autoregressive Moving Average Models", Jr.al. Amer. Statist. 

652-660, 1979. 

McKenzie, Ed. "General Exponential Smoothing and the Equivalent 
J. cf Forecasting, to appear, 193-a. 



McKenzie, Ed. "A Traditional Interpretation of the Forecasts cf Seasonally 
Differenced AP.IMA rrocesses* 1 , in Tine Aeries Aralvsist Thee rv end Pmc 
6, Proceedings cf ITS!.!, Toronto, August 19oj, Edited oy j.D .Ancersor., 



appear 193do 



Stem. G 



m 



series an mn 



( 4 i ^ 0 



ussier c: 



: crecastinr mi variate 






■ 1 1 "Kt j 



0 ) 0 ^ 
vo u 




Exhibit 1. US House soles, Jan 68-Dec 75. 

Solid line:Data; Dotted line:Current Level b Q l 
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Exhibit 2. US House sales: residual seasonal pattern 
after removal of Current Level b 0 l 
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